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Scaling tests in 0(a) -improved quenched lattice QCD* 
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We present a scaling investigation of renormalized correlation functions in 0(a)-improved quenched lattice 
QCD. As one observable the renormalized PC AC quark mass is considered, others are constructed such that 
they become the vector meson mass, and the pseudoscalar and vector meson decay constants in large volume. 
Presently, we remain in intermediate volume, (0.75 3 x 1.5) fm 4 , and study the approach to the continuum limit. 



1. Introduction 

Discretization errors of 0(a) in lattice QCD 
can be removed via a systematic approach based 
on the Symanzik improvement programme |Q, 
which adds appropriate higher-dimensional op- 
erators to action and fields of interest Ex- 
ploiting chiral symmetry restoration and certain 
current Ward identities on the lattice, a (mostly) 
non-perturbative 0(a)-improvement for action 
and quark currents as well as their renormal- 
ization has been achieved in the quenched case 
within this framework |j],|]J(|. 

Thus one is not only interested in the influ- 
ence of improvement on agreement of lattice data 
with experiment at fixed /3-values, but also in 
the quality of scaling and the size of its viola- 
tion. In this context it was reported || that at 
a ~ 0.1 fm the residual 0(a 2 ) lattice artifacts 
may be fairly large e.g. for JkTo (~ 10 %), while 
they are already very small for other quantities 
like m p /y/a (~ 2 %). Restricting to an interme- 
diate volume, we therefore examined the impact 
of 0(a) -improvement thoroughly and with high 
accuracy for different observables. 

2. Fermionic correlation functions 

Consider correlation functions in the 
Schrodinger functional (SF) Jl(| with all details 
found in H||. We use £(C) as xq = boundary 
(anti-)quark fields in O a = a 6 £ u v C(u) 75 ^C(v) 

and Q% = a 6 Eu,v C( u )7/cT"C(v), and axial (vec- 
tor) current A a ^(V^) and pseudoscalar (tensor) 
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density P a (T^ v ) to form the expectation values 



/a(zo) = -3 (A a (x)O a ) (la) 

Mzo) = -\(P a {x)O a ) (lb) 

h = -^(0' a O a ) (lc) 

fcv(xo) = -\{V£{x)Q%) (Id) 

M*o) = -UtU^QD- (le) 



Given the improvement coefficients ca and cy, 
the improved currents (Ai)^ = A^ + ac A d^P a 
and {Vi)l = V£ + ac v d„T^ lead to define the 
corresponding fermionic correlation functions 

f\{ x o) = /a(zo) +acAd fp(x ) (2a) 
k\/(xo) = fc v (xq )+ acy d k T (x ) , (2b) 

where the symmetrized lattice derivative acts as 
usual as d fl f(x) — [f(x + afi) — f(x — afi)]/2a. 

3. Observables under study and results 

Employing a mass-independent renormaliza- 
tion respecting 0(a)-improvement, the quantities 



(A a )« = Z A (l + 6 A am q )(Ai)° (3a) 
(Vk)£ - Z Y (l + 6 v am q )(y I )^ (3b) 
P£ = Z P (1 + b P am q )P a (3c) 

induce the renormalized correlation functions 

fJt(xo) = Z A (l + b A am q )f A (x ) (4a) 
ky-(x ) = Zv(l + bvam q )k v (x ) (4b) 



with Zx and bx, X=A,V,P, which are not func- 
tions of the quark mass am q = (1/k — 1/k c )/2. 
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Now we construct the following set of observables 
9o/a(zo) + ac A 9^/p(a;o) 



m(x ) 
mp S (xo) 
mv(x ) 
afpsM 



2f P (x ) 



dpfp(xo) 
fp( x o) 

dpfcy^o) 
ky(x ) 

V7T ' 

k^(x ) 



TO^r = mps(f ) 
m p = mv(f) 

/tt = Jps(-j) 



fp — /v ( 2 ) ' 



(5) 
(6) 
(7) 
(8) 
(9) 



The division by y^JT cancels the renormalizations 
of the boundary quark fields ensuring that / w ,/ p 
are scaling quantities, and the proportionality 
constants in eqs. (||),(^|) are such that these ra- 
tios turn, as T — > oo, into the familiar matrix ele- 
ments, which define the tt and p meson decay con- 
stants. The renormalized PCAC (current) quark 
mass in the SF scheme is obtained as M 



Za 
Zp 



topcac j "ipcAc = m(x ) 



(10) 



where (6a — 6p)am q may be neglected [pi. 

For the analysis we use c sw , cx and non- 
perturbatively determined in Hd,§| for (3 > 6.0, 
while 6a, &v and bp, as well as the SF specific 
improvement coefficients of the boundary coun- 
terterms ct and ct, are taken from 1-loop pertur- 
bation theory @J3JlT 



L/a 


/3 


K 


L/r 




8 


6.0 


0.13458 


1.490(6) 


2.004(9) 


10 


6.14 


0.13538 


1.486(7) 


1.946(14) 


12 


6.26 


0.13546 


1.495(7) 


2.050(16) 


16 


6.48 


0.13541 


1.468(8) 


1.991(15) 



Table 1. Simulation points for the LCP studied. 

The strategy was then to keep a finite physical 
volume and the quark mass fixed by prescribing 
the geometry T/L = 2, 'pion' mass m^L = 2.0 
and spatial lattice size L/ro = 1.49 using the re- 
cent results on the hadronic scale ro/a pi]] , see 
table [l]. This yields an intermediate volume of 
(0.75 3 x 1.5) fm 4 , moving on a line of constant 
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Figure 1. Scaling of the PCAC quark mass in in- 
termediate volume with SF boundary conditions. 
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Figure 2. Same as figure [I] but for the rho mass. 

physics (LCP) in parameter space with lattice 
resolutions ranging from 0.1 to 0.05 fm. As an 
important prerequisite for the reliability of the 
scaling test we verified by variation of the 1-loop 
values as simulation input that the dependence on 
ct, ct is small enough to be neglected. The leading 
scaling violations should therefore be 0(a 2 ). Any 
small mismatch with the renormalization condi- 
tions on m^L and L/tq was corrected by an es- 
timation of the corresponding slopes. Finally, we 
performed extrapolations to the continuum limit, 
assuming convergence with a rate cx a 2 . 



mro 


m p r 


Ur 


l/fp 


0.1092(45) 


1.860(20) 


0.704(13) 


0.826(22) 


2.1 % 


3.0 % 


4.8 % 


28 % 



Table 2. Continuum limits and their percentage 
deviations from (3 — 6.0 (a ~ 0.1 fm). 
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Figure 3. As figure [I] but for the ir decay con- 
stant. 
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Figure 4. As figure [1| but for the p decay con- 
stant. Diamonds refer to an analysis with cy = 0. 

The fits are displayed in figures |l| - |J, where 
the total error is always dominated by the uncer- 
tainties of the renormalization factors Zx- One 
observes the leading corrections to the contin- 
uum to be compatible with 0(a 2 ). Moreover, it 
can be inferred from table || that the difference of 
the continuum limits from the values at /3 = 6.0 
(a ~ 0.1 fm) is below 5 % in the improved theory. 
The only exception is the inverse 'rho' meson de- 
cay constant, whose slope is quite large; for that 
reason we discard the (3 = 6.0 point to extrapo- 
late to the continuum limit. 

4. Discussion and outlook 

Our numerical simulations of renormalized and 
improved correlation functions show an overall 
behaviour completely consistent with being linear 
in a 2 at (3 > 6.0 for all quantities under consider- 
ation. Changing a by a factor 2 gives very stable 



continuum extrapolations. 

But the residual 0(o 2 )-effect in l// p is still 
large at a ~ 0.1 fm (< 30 %). Here we find that 
artificially setting cy = results in data with 
an overall weaker dependence on a. However, as 
indicated in figure || for such a cy-value the func- 
tional form of the a-effects seems no longer com- 
patible with a 2 . Additionally, a chiral Ward iden- 
tity is badly violated for cy = at 0(a)-level M. 
Thus there is no choice of cy, which makes 0(a 2 )- 
effects small in both channels. This example 
clearly illustrates that even in the 0(a)-improved 
theory the remaining O(o 2 )-effects have to be as- 
sessed by varying the lattice spacing. 

This work is part of the ALPHA collaboration 
research programme. We thank DESY for allo- 
cating computer time to this project. 
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